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We develop a model for the evolution of economic entities 
within a geographical type of framework. On a square symme- 
try lattice made of three (economic) regions, firms, described 
by a scalar fitness, are allowed to move, adapt, merge or cre- 
ate spin-offs under predetermined rules, in a space and time 
dependent economic environment. We only consider here one 
timely variation of the "external economic field condition". 
For the firm fitness evolution we take into account a con- 
straint such that the disappearance of a firm modifies the 
fitness of nearest neighboring ones, as in Bak-Sneppen popu- 
lation fitness evolution model. The concentration of firms, the 
averaged fitness, the regional distribution of firms, and fitness 
for different time moments, the number of collapsed, merged 
and new firms as a function of time have been recorded and 
are discussed. Also the asymptotic values of the number of 
firms present in the three regions together with their average 
fitness, as well as the number of respective births and collapses 
in the three regions arc examined. It appears that a sort of 
critical selection pressure exists. A power law dependence, 
signature of self-critical organization is seen in the birth and 
collapse asymptotic values for a high selection pressure only. 
A lack of self-organization is also soon at region borders. 

PACS numbers: 789.65.Gh, 05.10.Ln, 89.75.-k, 07.05.Tp, 
05.65.-|-b 



I. INTRODUCTION 

The present contribution extends our earlier paper, 
[1,2]. World economic conditions evolved and are quite 
varied on different time and space scales. Basic questions 
are whether the consequences of political conditions have 
predictable effects, or not and whether annoying situa- 
tions can be avoided. The questions pertain to macroeco- 
nomic themes, not so often touched upon in econophysics 
(see exceptions in [3 5]). 

The case of the Berlin wall fall followed by Eastern 
Europe and Central Asia market openings to a so called 
liberal economy is a intriguing event. From an econo- 
physicist point of view, the event (s) can be considered 
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as an increase in "physical volume, or available space" 
as well as in a modification of the external fields [1]. In 
the latter reference wc have discussed whether one can 
describe, within a simple model, the concentration of en- 
terprises and their so called "fitness" as a function of 
time and space, under varying in time and space eco- 
nomic field conditions. We have observed a non-trivial 
behavior with cycle features; in our opinion more pro- 
nounced that those discussed by Aoki, Kalecki, Freeman 
and others [5-8]. The model however seemed to be too 
optimistic in not allowing enough bankruptcies, collapses 
or mere disappearance of firms in the process, in par- 
ticular with respect to the so called economic selection 
pressure [1,2]. Moreover wc neglected the fact that a firm 
disappearance could locally modify the probability of sur- 
vival of neighboring ones. Indeed an immense political 
worry concerns the effect of degradation of an economy 
through avalanche-like processes. 

In practice, adaptation to a dynamic environment con- 
tains a trade-off: it is more difficult to adapt in the next 
(time) step as shown e.g. also on financial time series by 
Yamasaki et al. [9]. It is also expected that the adap- 
tive behaviors of "populations" differ depending whether 
a Darwinian or Lamarckian scheme is implied. A large 
discussion exists whether such schemes hold or not, and 
how in economy, e.g. see [10-17] for references and a 
recent review. Alas, the definition of Lamarckism (or 
Darwinism) is not totally agreed upon in economy cir- 
cles. 

Therefore we present and discuss here below a more 
elaborate model than in Ausloos et al. [1,2], i.e. in order 
to describe changing basic economic conditions and im- 
plement more realistic Darwinian- like adaptation rules, 
as discussed in economy circles [10-17] - in space and 
time. In essence we introduce a Darwin-Bak-Sneppen- 
like constraint [18] in our previous macro-economy [1,2] 
microscopic physics-like model, i.e. the choice of the ini- 
tial firm or company depends on its relative fitness with 
respect to the external field as in [2]; the chosen firm is 
the less adapted one, instead of being chosen at random 
[1]. The fitness of neighboring ones evolve according to 
some a priori defined dynamics. 

The notion of economic external field is also important. 
Nelson and Winter [17] write that "at any time, firms in 
an industry can be viewed as operating with a set of 
techniques and decision rules (routines) keyed to condi- 
tions external to the firm" ... "and to various internal 
state conditions". Surely self-organization and external 
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(field) conditions exist at the same time, and are hard 
to separate [19,20]. However there arc also many dif- 
ferent conditions affecting many characteristic variables 
of a company or industry state. Thus searching for as 
simple as possible description and consequences of such 
external though often qualitatively only apparent con- 
straint, to be contrasted to self-organization is of inter- 
est. As usual in physics when intending extreme simpli- 
fication this field should be applied toward or coupled 
to some variable such that some function be optimized 
[21], like the free energy in thermodynamics. In macro 
and micro-economy the field can be imderstood as of po- 
litical origin, but also of more physical origin like the 
weather. In the latter case, see the weather influence on 
the energy prices [22] . As example of the former one can 
think about deregulation, or wage and benefit taxation 
influences upon localization and production or market 
approach by a company. 

Along this line, the most extreme simplification leads 
to represent a company by a variable having the same 
number of characteristic symmetry group elements as the 
field. Thus a scalar, called the fitness of a company is 
introduced [1,2]. One can consider that the fitness is 
the price of stocks and the external field can be in some 
sense the value of the S&P500 to which the price is com- 
pared. A buy /sell strategy will depend on the difference 
between the two values. A too large difference can sug- 
gest bankruptcy and disappearance of a company, or the 
'need" for a modification of the business plan. The notion 
of fitness of a company is related to such a book value. 
The fitness might also represent the number found on 
gross sales or benefit lines.... The field notion has to be 
subsequently adapted. 

Any such system will be evolutionary if it has a fourth 
ingredients: beside the set of degrees of freedom, each 
ranked by a quality fitness criterion, and a mechanism 
for introducing mutations to interject diversity into the 
productive process, the system must be governed by a 
selection process based on the fitness ranking [6,23,24] as 
in biology [25]. 

Thus we consider a system composed of firms located 
at the nodes of a locally square symmetry lattice of L^, Ly 
size dimensions. As in Ausloos et al. [1] the system is 
divided into three parts, k= I, II and III, by vertical 
boundaries at L^k, i-c lines along the y axis at Lxi^ Lx2- 
There is no periodic boundary conditions. At the be- 
ginning there arc N(t — 0) firms, located only in the 
first (I) region. One lattice site may be occupied by only 
one firm. Apart from its location on the lattice a firm 
i is characterized by its scalar fitness, /j, taken from a 
uniform distribution, e (0,1). There is no other factor 
distinguishing the firms, hence we may consider them to 
belong to only one group of industrial companies, i.e. 
representing a particular industrial branch. Initially the 
fitness values are attributed to firms in a random way. 

Any firm may change its position [20] to a vacant near- 
est neighboring site chosen at random and may either 
merge with another firm (with a given probability, - ac- 



cording to its business plan) or create a new one. In a 
region, a firm is under an external field F G (0,1), which 
characterizes the local economic conditions. The local 
political system may be more or less demanding,- more 
or less high taxes, salary rules, worker age (and other) 
conditions, this is taken into account through the pa- 
rameter, sel, - corresponding to selection in biological 
systems [25]. In our model it is also a scaling parameter 
for the field and the fitness, favoring or not the possible 
disappearance of a company unfit to its environment, the 
field [26-28]. 

The field F may (and will) differ from region to region 
and may also change in time, measured in Monte Carlo 
Steps (MCS). In the following we will only consider one 
drastic field change and keep stable (or static) all envi- 
ronmental (economic field) conditions for evaluating its 
effect on the population of entities evolution. 

As in Ausloos et al. [1] after some time the barrier 
separating the region I from regions II and III is open 
and the diffusing firms may enter the latter regions. At 
that time the external field also changes in the first region 
and some values are given to the fields in the regions II 
and III. 

Before any activity is taken for the firm (moving, merg- 
ing...) its condition is checked by comparing its fitness 
with the field F. If the two values differ strongly, the 
firm may collapse and be removed from the system (see 
the algorithm below). 

The Bak-Sneppen condition, i.e. a modification of the 
fitness of the eliminated firm nearest neighbors, is imple- 
mented through distributing new (uncorrclated in space 
and time) random fitness to these (at most 4) neigh- 
bors. This corresponds to e.g. simulating the effect a 
bankruptcy may have on other firms connected to the 
one which went down. Obviously the fitness neighbors 
can increase or decrease. Some firms may profit or regret 
the disappearance of a competitor in the local market. 

Firms will also be systematically moving on the lattice, 
looking for partners in view of merging or producing spin- 
offs. We will take into account the number of contacts 
in the company network (neighbors on a lattice in our 
case) for describing the firm evolution. It is empirically 
known that the replacement of some contacts is needed 
for innovation [26-29]. 

The (new) firm(s) fitness evolution rule is slightly dif- 
ferent from that in Ausloos et al. [1], but still takes into 
account the fitness of parents. 

The following quantities are recorded: concentration of 
the firms in the three regions and the averaged fitness in 
the three regions both as a function of time, the spatial 
distribution of firm concentration and fitness for different 
time moments, the number of collapsed, merged and new 
firms as a function of time, finally the asymptotic, i.e. 
when a stationary state is reached, values of the number 
of firms present in the three regions together with their 
average fitness, as well as the number of respective birth 
and collapses. 

In section 2, we outline the simulation technique used 
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for implementing the model. We present a few results in 
Sect. 3, and end with a conclusion in Sect. 4. 

II. MODEL 

The algorithm we implemented is the following one: 

1. from a number of firms N(t) at a given time t a 
firm j is randomly picked 

2. its survival probability is calculated from the con- 
dition 

Pj=exp{-sel\F-fj\) (1) 

and compared to a random number r taken from 
a uniform distribution £ (0,1). If r > pj, the firm 
collapses, i.e. it is removed from the system; the 
nearest neighbors of the eliminated firm receive new 
random numbers (taken from a uniform distribu- 
tion) for their fitness. The algorithm then goes back 
to 1, 

3. if r < , we try to move the firm to a nearest neigh- 
borhood (NN). A random number, n, is generated 
from a uniform distribution and if it is smaller than 
0.25 then we check whether the "Northern NN" is 
an empty site, and we move the firm there, if pos- 
sible; if ri is between 0.25 and 0.50 then we look 
into the Western NN, etc. If the displacement trial 
is not successful the algorithm goes back to 1 and 
search for a new firm; 

4. next we look for a partner in the nearest neighbor- 
hood of the displaced firm new position. If there is 

a firm in the NN, 

5. with a probability (5 =) 0.01 the neighbor, say i, 
merges with the old firm j, which changes its fitness 
to 

fj = 0.5(/i + + (0.5 - r2)\fi - fj\), (2) 

where r2 is a random number in the range (0,1). 
The firm i then disappears from the system; 

6. otherwise with a probability (1 — 6 =) 0.99 the firms 
i and j produce a new firm k, a spin — off. The k 
firm is randomly positioned in the Moore neighbor- 
hood of the j firm on an empty site if any exists. 
The procedure for finding an empty place is similar 
to the one used when looking for an empty site to 
move the firm, except that on a square lattice the 
Moore neighborhood consists of 8 sites - NW, N, 
NE, W, E, SW, S and BE; whence if the random 
number ra is smaller than 0.125 we check the site 
NW, if it is larger than 0.125 but smaller than 0.250 
we check the N site, etc. The new firm receives its 



fitness depending on that of both parents (which 
remain in the system) as 

fk = 0.5(/i + + (0.5 - r4)|/i - fj\), (3) 

where and r4 are random numbers in the range 
(0,1). 

7. When N{t) firms were picked through (1), one MCS 

is said completed. 

Typical values taken for the discussed simulations be- 
low were : L^i = 50, = 100, Ly = 201, c(0) ~ 0.8, 
corresponding to a number of firms A^(0) = 8040. 

The parameters of the model are: sel, selection pres- 
sure; tchangej time after which the barrier at L^i is open 
and the field changed in I; the values of the fields. At 
the beginning the field in region I has the value Fj = 
0.5, which is "optimal"; the values in the regions II and 
III are irrelevant before the opening of the L^i barrier, 
since there are no firms there. After the change, at time 
tchange = 100 MCS, the valucs were arbitrarily taken to 
be: Fj = 0.3, F// = 0.5, Fm = 0.6, which means that 
the conditions in the region I deteriorated and the best 
situation is in region II. The conditions in region III are 
also imposed to be better than those in region I. These 
values remain unchanged till the end of the simulations. 

Although the presented curves were obtained each 
from a single simulation, we have checked that a different 
initial distribution of the firms (for the same initial con- 
centration, external field and selection) leads to a very 
similar, even quantitatively, set of situations. We are not 
interested here in the transient regime before the first 
100 MCS. The simulations were carried out till a station- 
ary state was reached, when each investigated quantity 
mildly oscillated around a rather stable average value. 
This happened quite soon (several hundreds of MCS) for 
low selection, but we had to run till 20000 MCS for the 
highest selection seA = 1.7. For a sel as "high " as 1.8, 
we found that no firm survived after the initial 800 MCS 
in the whole system. 

III. RESULTS AND COMMENTS 

There are several ways of presenting pertinent results. 
In the following we stress cases demonstrating the perti- 
nence of the model in view of economic qualitative obser- 
vations. The scales are chosen in order to allow for better 
visually comparing cases and emphasizing differences in 
behavior. 

Recall that from time tchange = 100 MCS, the field 
values are kept to be Fj = 0.3, Fjj = 0.5, and Fm = 
0.6 respectively. 
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A. Time dependent and regional values 

In Fig. 1 (a-b) the number of new firms appearing in 
each region is shown as a function of (MC) time, (a) in 
the case of a rather low sel, i.e. = 0.7 and in the case (b) 
of a rather high sel = 1.7. In Fig. 2 (a-b) the number 
of eUminatcd firms versus time, for the above two values 
of sel are shown. It is worth pointing here, the differ- 
ence with respect to the case examined in ref. [1] . In the 
latter paper the long time stability (asymptotic value) 
was obtained through a sharp increase followed by an 
exponential decay toward a zero death (and also birth) 
evolution. In the present case, the number of born or 
eliminated firms increases smoothly and remains stable 
at a finite value. Recall that there arc 8000 firms ("possi- 
ble parents") at t=0. Observe the x-axis scale indicating 
the effect of selection pressure on first time invasion of 
regions II and III. 

A minimum in the birth and death number is always 
observed at short time (after 100 MCS) in region I. No- 
tice from Fig. 1(b) and Fig. 2(b) that a strong sel might 
quickly kill any evolution process : zero-birth and zero 
death. It is also clear from comparing Fig. 1 and Fig. 
2 that there is a rapid increase in concentration in the 
newly opened regions as soon as the barrier at L^i is re- 
moved. The death and birth processes are always quan- 
titatively equivalent, indicating that a detailed balance 
equilibrium is intrinsically met in the model. 

The number of firms present in the three regions vs. 
time (in MCS) is shown for several sel values, from 0.4 
till 1.7 in Fig. 3 (a-d). At very low sel the number of firms 
rather quickly reaches a rather high saturation value in 
each region, in fact approximately equal to the initial 
value in region I. At high sel, Fig. 3 (c-d), there is a 
marked drop in region I, while regions II and III are in- 
vaded much later than in Fig.3 (a); the saturation level 
has a lower value than the initial starting one in region 
I. The asymptotic number of firm values is seen to de- 
crease with sel. The strongly marked oscillations might 
be considered as a signature of cycles [5,8,24]. At high 
sel values the effects are much enhanced. Also, regions 
II and III are invaded very slowly and late in time. Yet 
the fitness optimal value in each region is easily reached 
and remain stable in time (see Fig. 4 below). 

In Fig. 4(a-d) we show the average fitness / in the 
three regions vs. MCS time for several values of the se- 
lection pressure, ranging from 0.4 to 1.7. Recall that 
the values of the external field are respectively 0.3, 0.5 
and 0.6, and the wall at L^i falls at t=100. For a low 
sel (=0.4) small fluctuations are seen around the aver- 
age (asymptotic) values which are close to 0.5 in all re- 
gions. Since the selection pressure is low, the firms are 
not strongly forced to being in perfect matching to the 
field value, whence large deviations from the regional av- 
erage fitness can be enforced. 

For higher sel, data is quite scattered, and indicate a 
strongly hard learning (adaptation) process, but the av- 



erage fitness in region I is near 0.3, and that in region 
III close to 0.6 as expected. For high sel large fluctua- 
tions in the average / value are visible as in Fig. 4(b), 
at rather short time. At very large sel the beginning of 
the process is rather noisy, but the average fitness stabi- 
lizes at a value expected from the imposed external field. 
We adopted scales for the figures in order to show qual- 
itatively the decrease in the variance of the fitness dis- 
tribution with time, for various selection pressures. The 
adaptation is clearly faster in region I and II, and the 
variance much reduced as soon as the expected average 
values are obtained. 

Considering the physical diffusion process, we recorded 
the maximum distance reached by the right most firm as 
a function of time, and present it on a log-log scale for 
different sel values, from 0.3 till 1.8, in Fig. 5. Recall 
that due to the lattice topology, critical values are to be 
found at x= 50, 100 and 150. Studying Fig. 5, starting 
from high sel values, it appears that a plateau occurs (for 
sel =1.8 and 1.5) near _Lj,2=100, also somewhat seen for 
sel = 1.1. The ffow into region II and III seems more 
continuous at small sel values. The small tails sticking 
upwards (above x= 50) for t > 1000 witness unsuccessful 
attempts to colonize regions II and III, suggesting that 
a se/-dcpendent nucleation-growth process exists in the 
model. Moreover, for high sel values it appears that the 
region near the border at Lxi= 50 is sometimes depleted., 
indicating a gain that high sel values might destroy the 
population at all. From Fig. 5 (and others not shown 
here) we can calculate the time necessary for the first 
firm to reach the region right boundary, i.e. respectively 
-^x2=100 and L^s = 150, as a function of the selection 
pressure (Fig. 6). The curves look like high order power 
laws (with exponents equal ca. 4 and 5 respectively), but 
more complicated functions might be tried. 

B. Asymptotic values 

It is of interest to observe beyond the transient regimes 
the equilibrium values. In Fig. 7 (a) the normalized 
asymptotic values of the born firms in the three regions 
are shown as a function of sel. The "normalization" is 
made with respect to the total asymptotic number of 
firms in the region shown in Fig. 7(b). "Asymptotic" 
meaning here that we have averaged the last 500 entries 
for 10 000 MCS runs. The same type of data is obtained 
quanfAtatively for the number of disappearing firms (not 
shown). A maximum occurs near sel = 0.8 with different 
analytic behaviors on both sides of this value. The in- 
crease goc^s lik(^ a stretched exponential as for diffusion- 
controlled growing entities [30] at low sel (< 0.8), i.e. 
sel^'^^exp{—0.5/ y^sel)), but the decay at high sel has 
been found to be markedly a power law sel~°''^'^, like for 
critical systems [31]. Recall that each region may contain 
at most 10 050 entities, a number to be compared to the 
number of asymptotically existing firms, as displayed on 
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Fig. 7 (b). For completeness, the asymptotic fitness in 
the three regions is shown vs. sel in Fig. 8. This indi- 
cates that only for large sel (> 0.8) has the asymptotic 
average fitness almost reached the expected one from the 
field conditions. 



C. Regional behaviors 

We show the short time dependence {t <800 MCS) of 
the (average over a "vertical column" ) firm concentration 
vs. the column number "a;" for rather strong (1.3) and 
very strong (1.7) sel values in Fig. 9. A drop in the over- 
all concentration is well seen in region I at short time, but 
the concentration recovers later on, as seen from previ- 
ous paragraphs and figures. For sel =1.8 the selection 
pressure is too high to maintain any entity in region I 
after a few MCS. 

The dependence of the firm (vertical) concentrations 
on position " x" at different (longer) times is shown in 
Fig. 10 (a-c) for sel values just above the critical one, 
i.e. ~ 0.8. 

Interesting and unexpected features occur. At very low 
sel (not displayed) a saturation level is reached right from 
the beginning of the simulation, indicating "easy adapta- 
tion" everywhere. At high sel the third region is invaded 
ca. after 5000 MCS, while region II is invaded after 800 
MCS. We emphasize the occurrence of large well marked 
"oscillations" in the local concentrations. Dips do occur, 
see Fig. 10 near the region borders, - likely due to the 
difficulty for firms to reach the most appropriate fitness 
together with their neighbors, due to the interregion field 
gradient. 

Finally the (vertically averaged) fitness vs. position 
along the a:-axis taken at three time steps and for two 
sel values on both sides of the critical one are shown 
in Fig. 11. Notice the different time(s) chosen for such 
snapshots. For medium sel (=0.7), the scatter around 
the expected theoretical fitness imposed by the external 
field is not large but the average values nevertheless do 
not closely correspond to the expected ones from the im- 
posed field, in all three regions, - the more so on the 
average , and asymptotically as seen above. At high sel 
(=1.7), the data is more scattered but the expected val- 
ues are better recognized as fulfilling the external field 
condition. 



IV. CONCLUSION 

The dynamics of an economic firm population has been 
considered through a model considering entities charac- 
terized by a scalar number, diffusing on a regular lattice, 
merging, collapsing or creating spin-offs. The economic 
environment is described by a static regionally dependent 
field. The adaptation and evolution of the firms to the 
field condition have been mimicked by a Darwinian-like 



selection rule but with a fitness evolution equation for 
the newly appearing firms taking into account that of the 
parents. The firm fitness distribution is shown to have 
some best adaptation difficulty, as discussed by Yamasaki 
et al. [9]. We have found an unexpected and interesting 
feature: the local changes of the environment is leading to 
sharp variations, almost discontinuous ones, in the fitness 
and concentrations, in particular when the field gradient 
is strong. A lack of self-organization is thus seen at re- 
gion borders. This is due to the fact that for an entity 
attempting an adaptation, the learning process is quite 
hill climbing as in NK models [32]. Indeed the evolution 
of an economy, in which the functioning of companies 
is interdependent and depends on external conditions, 
through natural selection and somewhat random muta- 
tion is similar to bioevolution on NK fitness landscapes as 
described by Kaufmann [32]. In these evolutionary eco- 
nomics models [28] , economic agents randomly search for 
new technological design by trial-and-error and run the 
risk of ending up in sub-optimal solutions due to inter- 
dependencies between the elements in a complex system. 
As argued by Prenken [28] these models of random search 
are legitimate for reasons of modelling simplicity, but re- 
main limited as these models ignore the fact that agents 
can apply heuristics." We totally agree. Indeed there 
is no mimicking nor endogeneous search for optimisation 
as in [26-28] in the present model, the algorithm being 
only geared toward reducing the distance between the 
company fitness and the external field value. 

One might take into account in a more proper way the 
number of business contacts, the evolving firm neighbors 
: the development of a market or a business plan, or a 
macroeconomy depends on the available number of con- 
tacts (bonds between nodes on a network) whence to the 
creation and destruction of contacts. A large stock of re- 
lational capital (business contacts) usually increases the 
sold output of a representative firm. Other lattices or 
network structures should be usefully studied. 

It would be of interest to consider specific (historical) 
cases, and connect the MC simulations time scales to 
real cases. The time scale of changing environments is 
an available parameter which could indicate how robust 
the model can be with respect to realistic adaptation and 
evolution. Some flexibility for defining time scales exist 
in the (i) b and (ii) sel values: (i) The value of the spin- 
off creation probability might be too large for previous 
economies. It has been really shown in the line of the 
ACP model that for the best fitted company a small b 
value leads to chaotic behavior while a large b leads to a 
monopoly like situation [24]. (ii) A sel (space and time 
dependent) value could be found by examining macroe- 
conomy data. A sort of critical selection pressure exists 
separating different birth-death regimes. This might be 
used as a time and field scale definition when counting 
bankruptcies in countries, including per industry type. 

Contrary to the Bak-Sneppen model [18], the average 
fitness(es) or concentration(s) do not seem to appear as 
power laws (like those characterizing avalanches), except 
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in the high selection pressure regime and for long (asymp- 
totic) times. The power law dependence, which is a usual 
signature of self-critical organization is only seen for a 
high selection pressure in the asymptotic behavior of the 
birth-death number. Nevertheless from the set of results, 
as presented here above, we observe that in such a Dar- 
winian evolving economic world: (i) there are relatively 
well marked effects due to the "selection pressure"; (ii) 
temporal field changes can imply a stable density distri- 
bution, as in Ausloos et al. [1]; (iii) a diffusion process 
together with a business plan, and the selection pressure 
smoothly leads to asymptotic equilibrium states with re- 
spect to birth and death processes; (iv) the fitness does 
not always reach the a priori externally imposed field 
value; (v) the role of the gradient (economic) field at bor- 
ders might indicate either complex oscillations or chaotic 
behaviors in such regions. 

We are aware that further improvements are needed. 
We are drastically caricaturing macro and micro economy 
field conditions, as well as the description of the "inter- 
nal" interactions sequence(s). Surely a company, or a set 
of industries, should not be described by one scalar num- 
ber fi, but rather a vector (or matrix) model coupled 
to a (so called external) vector (or matrix) field should 
be examined. Moreover the birth and death process de- 
scription through merging and spin off's could also be 
stochastic or evolutive. The spatial distribution in the 
distinct regions might be also on interest [29] as in biol- 
ogy [33]. 
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Fig. 7 - (a) Normalized asymptotic values values of the 
born firms in the three regions as a function of various sel; 
"Normalization" is with respect to (b) total asymptotic 
number of firms in the respective region 

Fig. 8. - Asymptotic fitness in regions I, II and III vs. 
sel ; recall that the field values are imposed to be 0.3, 05. 
and 0.6 respectively 

Fig. 9 - Short time dependence of the firm concentra- 
tion in regions I, II, III vs. spatial coordinate ".t" for sel 
values (a) = 1.3, (b) = 1.7 and at different time moments 
- 10,20,50 and 800 MCS. 

Fig. 10 - Dependence of the firm concentrations on 
position "x" at different times for sel = (d) 1.1; (e) 1.5; 
(f) 1.7 

Fig. 11 - Vertically averaged fitness vs. position along 
the a;-axis at three time steps and for two sel values, i.e. 
(a) 0.7 and (b) 1.7 



FIGURE CAPTIONS 



Fig. 1 - Number of new firms born in region I, II, III ( 
-|-, X, dotted square, respectively) as a function of time, 
(a) sel = 0.7 ; (b) sel = 1.7 

Fig. 2 - Number of eliminated firms in the three re- 
gions versus time, for (a) sel = 0.7 ; (b) sel = 1.7 

Fig. 3 - Number of firms in the three regions vs. time 
(in MCS) for (a) sel = 0.4; (b) sel = 0.9; (c) sel = 1.3; 
(d) sel = 1.7 

Fig. 4 - Average fitness in the three regions I, II, III 
( -|-, X, dotted square, respectively) vs. MCS time for 
several values of the selection pressure (a) = 0.4; (b) = 
0.9;(c) = 1.3; (d) = 1.7 

Fig. 5 - Maximum distance reached by the right most 
firm as a function of time, on a log-log scale for different 
sel values 



Fig. 6 Shortest time necessary for a firm to reach the 
region II (100) or III (150) right boundary, i.e. Lxi=100, 
-^a;2=150, as a function of the selection pressure ; curves 
indicate the best power law fit, i.e. ca. 4 (top) and 5 
(bottom) respectively 
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